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Abstract 



We calculate the cross section of the inverse beta decay process, + n ^ p + e, in a, 
, background magnetic field which is much smaller than rrip/e. Using exact solutions of the 

^ I Dirac equation in a constant magnetic field, we find the cross section for arbitrary polarization 

CO . of the initial neutrons. The cross section depends on the direction of the incident neutrino 

I even when the initial neutron is assumed to be at rest and has no net polarization. Possible 

• implications of the result are discussed. 

o ■ 
<^ 

: 1 Introduction 

Q-i! The interactions of elementary particles show novel features when they occur in non-trivial 

Q^l backgrounds. Study of particle propagation in matter has proved pivotal in the understanding 

of the solar neutrino problem. Similar studies of particle processes in background magnetic 
fields are also important since stellar objects like neutron stars are expected to possess very high 
magnetic fields, of the order of 10^^ G or higher. Analysis of these processes might be crucial 
^ I for obtaining a proper understanding of the properties of these stars. 

' In this paper, we calculate the cross section of the inverse beta-decay process in a magnetic 

field. We consider the possibility that the neutrons may be totally or partially polarized in the 
magnetic field, and find the cross section as a function of this polarization. The neutrinos are 
assumed to be strictly standard model neutrinos, without any mass and consequent properties. 
The presence of the magnetic field breaks the isotropy of the background, and a careful cal- 
culation in this background reveals a dependence of the cross section on the incident neutrino 
direction with respect to the magnetic field. 

Considerable work has been done on the magnetic field dependence of the URCA processes 
which have neutrinos in their final states ^ [21 El |^ |5] . An angular dependence obtained in 
the differential cross section of these reactions imply that in a star with high magnetic field, 
neutrinos are created asymmetrically with respect to the magnetic field direction. The process 
that we consider, on the other hand, have neutrinos in the initial state. So this process influences 
the neutrino opacity in a star. 
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Some calculations of this process exist in the literature. Roulet [H], as well as Lai and Qian 
[2] performed the calculation by assuming that the magnetic field effects enter only through the 
phase space integrals, whereas the matrix element remains unaffected. Gvozdev and Ognev [2j 
considered the final electron to be exclusively in the lowest Landau level. Arras and Lai |1] 
calculated only the angular asymmetry, and only to the first order in the background magnetic 
field. The earlier calculations of the present authors [HI El did not take neutron polarization into 
account. In this paper, we consider the problem in full detail — i.e., we calculate the matrix 
element using spinor solutions of the electron in a magnetic field, take all possible final Landau 
levels into account, include the possiblity of neutron polarization, and perform the calculations 
to all orders in the background field in the 4-fermi interaction theory. 

The paper is organized as follows. In Sec|2 we provide some background for the calculation. 
Most of this section contains no original material, but we provide it for the sake of completeness, 
as well as for setting up the notation that will be used in the later sections. In Sec. 01 we define 
the fermion field operator and show how it acts on the states in the presence of a magnetic field. 
Sec. jl] contains the calculation of the cross section for a monochromatic neutrino beam, which 
contains the main results of our paper. In Sec. we discuss the realistic case where the initial 
neutrino beam has a finite energy spread. Sec. El contains our discussions and conclusions. 

2 Solutions of the Dirac equation in a uniform magnetic field 

For a particle of charge eQ, the Dirac equation in presence of a magnetic field is given by 



where (j) ^'^id x are 2-component objects. We use the Pauli-Dirac representation of the Dirac 
matrices, in which 



where each block represents a 2 x 2 matrix, and a are the Pauli matrices. With this notation, 
we can write Eq. 1)2. 1|) as 











(2.5) 



Eliminating X; we obtain 




(2.6) 
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We will work with a constant magnetic field B. Without loss of generality, it can be taken along 
the z-direction. The vector potential can be chosen in many equivalent ways. We take 



Ao = Ay = A^ = 0, 
With this choice, Eq. (|2.6j) reduces to the form 



A, 



-yB. 



m 



d 



V" + {eQB)Y - eQB{2iy— + a,) 



(2.7) 



(2. 



Noticing that the co-ordinates x and z do not appear in the equation except through the deriva- 
tives, we can write the solutions as 



(2.9) 



where f{y) is a 2-component matrix which depends only on the y-coordinate, and possibly some 
momentum components, as we will see shortly. We have also introduced the notation X for the 
spatial co-ordinates (in order to distinguish it from x, which is one of the components of X), and 

for the vector X with its y-component set equal to zero. In other words, p-X>^= p^x + pzZ, 
where px and pz denote the eigenvalues of momentum in the x and z directions.^ 

There will be two independent solutions for /(y), which can be taken, without any loss of 
generality, to be the eigenstates of with eigenvalues s = ±1. This means that we choose the 
two independent solutions in the form 



Uiy) 



f-iv) 



F-{y) 



Since azFg = sFg, the differential equations satisfied by Fg is 



dy'^ 



{eQBy + pxYFs + {E' 



m 



pi + eQBs)Fs = , 



(2.10) 



(2.11) 



which is obtained from Eq. (|2.<S|) . The solution is obtained by using the dimensionless variable 

Px 



^ = ^e\Q\B{y 
which transforms Eq. (|2.1ip to the form 



eQB 



(2.12) 



where 



E' 



0, 



m 



p1 + eQBs 



(2.13) 

(2.14) 
2v + 1 for 

(2.15) 

^It is to be understood that whenever we write the spatial component of any vector with a lettered subscript, 
it would imply the corresponding contravariant component of the relevant 4-vector. 



e\Q\B 

This is a special form of Hermite's equation, and the solutions exist provided 
1/ = 0, 1, 2, • • •. This provides the energy eigenvalues 

E"^ =m^ +pI + {2v + l)e\Q\B - eQBs , 
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and the solutions for are 



(2.16) 



where H^, are Hermite polynomials of order z/, and N^, are normalizations which we take to be 

1/2 



(2.17) 



We stress that the choice of normalization can be arbitrarily made, as will be clarified later. 
With our choice, the functions 1^ satisfy the completeness relation 



V 

where is obtained by replacing y by in Eq. (|2.12l) . 

So far, Q was arbitrary. We now specialize to the case of electrons, for which Q 
solutions are then conveniently classified by the energy eigenvalues 

E\ = n? +v\^ 2neB , 



(2.18) 



-1. The 



(2.19) 



which is the relativistic form of Landau energy levels. The solutions are two fold degenerate in 
general: for s=l, z^ = n— 1 and for s = —1, v = n. In the case of n = 0, only the second 
solution is available since v cannot be negative. The solutions can have positive or negative 
energies. We will denote the positive square root of the right side by Representing the 
solution corresponding to this n-th Landau level by a superscript n, we can then write for the 
positive energy solutions, 



f^-\y) 



(2.20) 



For n = 0, the solution /+ does not exist. We will consistently incorporate this fact by defining 

l-i{y) = 0, (2.21) 

in addition to the definition of /„ in Eq. (|2.1(i|) for non-negative integers n. 

The solutions in Eq. ()2.2U|) determine the upper components of the spinors through Eq. 
1)2. 9|1 . The lower components, denoted by x earlier, can be solved using Eq. (|2.5|) . and finally 
the positive energy solutions of the Dirac equation can be written as 



e-^P-^^C/,(y,n,a), 
where denotes the space-time coordinate. And Us are given by 

In-liO 





(2.22) 



/ 



\ 



Pz 



En + m 
V2rieB 



IniO 



( 



\ En + m J 





V2neB 



En + m 

Pz 



\ En + m 



(2.23) 
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A similar procedure can be adopted for negative energy spinors which have energy eigenvalues 
E = —En- In this case, it is easier to start with the two lower components first and then find 
the upper components from Eq. (|2.4|) . The solutions are 



where 



Pz 



En + m 
V2rieB 



V 



En + m 




IniO 



\ 



(2.24) 



( ^ ,~\ 

hin + m 

Pz 



J 



En + m 



In{i) 



V 







(2.25) 



where ^ is obtained from ^ by changing the sign of the p^-term. 

For future use, we note down a few identities involving the spinors which can be obtained 
by direct substitutions of the solutions obtained above. The spin sum for the ?7-spinors is 

Pu{y,y*,n,'^) = '^Us{y,n,'^)Us{y^,n,^,) 

s 

1 



2{En + m) 



\m{l + az) +i)\\ -|^||75|/„_i(^)/„_i(^^ 

+ [ m(l- a,)+i)\\ +^||75}/n(0^(e0 
\/2^(7i - i72)/n(04-i(CJ 
V2neS(7i + f72)/„-i(0^(C* 



(2.26) 



where we have introduced the following notations for any object a carrying a Lorentz index: 

--^ - (ao,0,0,a^) , 



ay 

sji* = (0^,0,0,00) . 

Similarly, the spin sum for the l^-spinors can also be calculated, and we obtain 

Pv{y,y*,n,^,) = J2ys{y,n,i^)Vs{y 

s 

1 



(2.27) 



, n, 



2iEn + m) 



+ 



{-m(l + a,)+|^ll - 1^1175} /n-i(e)^-i(e^ 

{-m(l - a,)+^i\ +j>l\l5} IniOlni^*) 
+^/2^(7i - ^72)^(0^-1 (6.) 

+\/2^(7i + i72)i„_i(e)/n(e^) 



(2.28) 
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3 The fermion field operator 



Since we have found the solutions to the Dirac equation, we can now use them to construct the 
fermion field operator in the second quantized version. For this, we write 



s=±n=0 



dpx dpz 
D 



fsin,ji)e-'P-''ms{y,n,ii) + p,in,ii)e'P-''Wsiy,n,ji)] . (3.1) 



Here, fsin^p^) is the annihilation operator for the fermion, and fl{n,fQ is the creation operator 
for the antifermion in the n-th Landau level with given values of px and pz ■ The creation and 
annihilation operators satisfy the anticommutation relations 



fs{n,p,),fl,{n,p^) = 6ss'Snn'S{Px -p'x)S{pz -p'z) , 



(3.2) 



and a similar one with the operators / and P, all other anticommutators being zero. The 
quantity D appearing in Eq. (|3.1|) depends on the normalization of the spinor solutions, and 
this is why of the spinors could have been chosen arbitrarily, as remarked after Eq. (|2.17|) . Once 
we have chosen the spinor normalization, the factor D appearing in Eq. (|3.H) is however fixed, 
and it can be determined from the equal time anticommutation relation 



(3.3) 



Plugging in the expression given in Eq. (|3.1j) to the left side of this equation and using the 
anticommutation relations of Eq. (|3.2() , we obtain 



dpx dpz 

L>2 



(^e-^M--*)e-'P^(^-^*^Usiy,n,j^)UUy.,n,ii) 
+e^P^(--*)e^P^(^-^*)Vs{y,n,fi)vHy.,n,ji)) . (3.4) 



Changing the signs of the dummy integration variables px and pz in the second term, we can 
rewrite it as 



EE 



dpx dpz 

L»2 



-ipx{x-Xi,) -ipz{z-Zi,) 



(Usiy,n,p^)Ul{y^,n,p,) 



+Vs{y, n, -p,)V}{yi,,n, -p,) 



(3.5) 



Using now the solutions for the U and the V spinors from Eqs. p. 23)1 and (|2.25)) . it is straight 
forward to verify that 

E {Us{y,n,p^)Ul{y^,n,p^) +Vs{y,n,-p,)V}{y^,n,-p^ 

s 

p1 + 2neB\ 



1 + 



(En + m) 







X diag i„_i(e)/n-i(c.),/n(0^(^*),^-i(0^-i(e^),^(e)/n(e.) (3.6) 



where 'diag' indicates a diagonal matrix with the specified entries, and ^ and involve the same 
value Px- At this stage, we can perform the sum over n in Eq. (|3.5|) using the completeness 
relation of Eq. ()2.18|) . which gives the 5-function of the y-coordinate that should appear in 
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the anticommutator. Finally, performing the integrations over px and pz, we can recover the 
5-functions for the other two coordinates as well, provided 

1 1 . ^ 

(3.7) 



En + m L>2 (27r)2 ' 

using the expression for the energy eigenvalues from Eq. 1)2. 19() to rewrite the prefactor appearing 
on the right side of Eq. 1)3. 6|) . Putting the solution for D, we can rewrite Eq. ()3.1|1 as 



oo „ 

^(^) = E E / 

s=±n=0 



dpx dpz I En + m 
2-K V 2Er, 



X [/,(n,;^)e-^f-^^C/,(2/,n,a)+/I(n,|^)e^P-^vy,(y,n,;?^)] . (3.8) 

The one-fermion states are defined as 

|n,:?^) = C7/t(n,g)|0) . (3.9) 

The normalization constant C is determined by the condition that the one-particle states should 

be orthonormal. For this, we need to define the theory in a finite but large region whose 
dimensions are L^, Ly and along the three spatial axes. This gives 

C = -^. (3.10) 

Then 



^l^u{X)\n,p,) = W ^" + 7 e-^P-''ms{y,n,p,)\0) , (3.11) 
where ipjj denotes the term in Eq. (|3.8() that contains the [/-spinors. Similarly, 



(n,a|v^f,(X) =^-^±^e^f^C/.(y,n,a)(0| . (3.12) 

4 Inverse beta-decay 

In this section, we calculate the cross section for the inverse beta-decay process Ve + n ^ p + e~ 
in a background magnetic field. The magnetic field might provide a net polarization of the 
neutrons, which we take into account. However, the magnitude of the field is assumed to be 
much smaller than m^/e or m^/e, so we ignore its effects on the proton and neutron spinors. 
The electron spinors, on the other hand, are the ones appropriate for the Landau levels. Thus, 
we can write the process as 

Ve{k) + n{P)^p{P') + e{ji,n'). (4.1) 
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4.1 The S-matrix element 

The interaction Lagrangian for this process is 

^int = \/2G/3 V(e)7^-^V'(i.e)J [V'(p)7m(5V -5A75)V'(n)J , (4-2) 

where L = ^(1—75) and Gp = Gp cos 9c, Oc being the Cabibbo angle. In first order perturbation, 
the S'-matrix element between the final and the initial states of the process in Eq. (|4.1|1 is 
therefore given by 



S 



fi 



V2Gp J 



(fx (e(^,n) V'{e)7''^V'(i.e) ^e{k) 



X ^p(P') V'(p)7A.(5y -5A75)V'(n) n{P)) . (4.3) 

For the hadronic part, we should use the usual solutions of the Dirac field which are normalized 
within a box of volume V , and this gives 

f,i{P'-P)X 



P{P ) iJ(p)lti{9v - 9Al5)ip{n) n{P) 



'^{p)iP'hf^i9V - 9Alb)U{n){P) ,(4-4) 



using the notations £ = Pq and £' = Pq. For the leptonic part, we need to take into account the 
magnetic spinors for the electron. Using Eq. 1)3. 12() . we obtain 



m 



Putting these back into Eq. 1)4. 3|) and performing the integrations over all co-ordinates except 
y, we obtain 



Sfi = {2^f^{P + k-P'-p') 



En' + m 



1/2 



(4.6) 



_2ujV 2£V 2£'V2En'L^L^_ 
Here, 6^ implies, in accordance with the notation introduced earlier, the (5-function for all space- 
time co-ordinates except y. Contrary to the field-free case, we do not get 4-momentum conser- 
vation because the y-component of momentum is not a good quantum number in this problem. 
The quantity is the Feynman amplitude, given by 



= y2G^[n(p)(P')7/.(5y - 5A75)n(„)(P)J J dy e''^yy[U(^,^{y,n' ,i()j''Lu^,^^){k)\ , (4.7) 
using the shorthand 

Qy = Py + h-Py- (4-8) 
The transition rate in a large time T is given by \Sfi\^/T. From Eq. KEii . using the usual 



rules like 



we obtain 



5{£ + uj - £' - En' 
6{P, + k,-P^-p'^ 
5{P, + k,-Pi-p'^ 

^ 3 r3/ 



' = Ls{£ + uj-£' -En'), 
Ztt 

' = ^^(p-+k,-p:,-p',). 



SMVT = -{27:r&%P + k-P' -p' 



En' + m 



16 



V^Lo££'En' 



(4.9) 
(4.10) 
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4.2 The scattering cross section 



Using unit flux 1 jV for the incident particle as usual, we can write the differential cross section 
as 



da = V 



dp, 



where dp, the differential phase space for final particles, is given in our case by 

V 



Therefore 



da = V 



dp 



P' . 



(4.11) 



(4.12) 



dp'x dp'z 



V 



d^P' 



' ^Hp + k-p'-p')^"^'^"^ 



647r2 



The square of the matrix element is 



UjSS'En' 



'^dp'^dp'zd'P'. 



(4.13) 



(4.14) 



where Hfj_y is the hadronic part and H.^'^ the leptonic part, whose calculation we outline now. 

For the hadronic part, we can use the usual Dirac spinors because of our assumption that 
the magnetic field is much smaller than m^/e. We will work in the rest frame of the neutron. 
Due to the presence of the background magnetic field, the neutrons may be totally or partially 
polarized. We define the quantity 



S = 



(4.15) 



where A^^^ denote the number of neutrons parallel and antiparallel to the magnetic field. Then 



i(i + 5)ijw + i(i-s)/f^;), 



(4.16) 



where H^fui^ denotes the contribution calculated with spin-up and spin-down neutrons respec- 
tively. Either of these contributions can be calculated by using the spin projection operator, 
which is ^(1 ± 7573) for up and down spins. A straight forward calculation then yields 



H^, = 2{gl + g\){P^,Pl + P^P'^ - g^^uP ■ P' 



+2{gy - gAjmniripg^^ + ^igv9A£fMuXpP P''' 
- S 4.gvgAmn{Plgzy + Plg^ix - P^g^v) + 2ie^,„2.aR'^ 

where we have introduced the shorthand 

R"" = {9v + gDmnP"" - {gl - g\)mpP'' . 



(4.17) 



(4.18) 
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We have omitted some terms in the expression for H^y that involve spatial components of the 
neutron momentum, with the anticipation that we will perform the calculation in the neutron 
rest frame. 

In the leptonic part , we should use the magnetic spinors given in Sec. El This gives 



dy / dy, e^<iy^y*-v) Tr y., n 



(4.19) 



where P\j denotes the spinor sum for the electrons, given in Eq. 1)2.26(1 . We now have to perform 
the integrations over y and y^. Each of these variables should be integrated in the range —■^Ly to 
+^Ly. However, since we will take the infinite volume limit at the end as usual, we let L^^ — > oo 
and use the result ^ 



This gives 



where 



1 



eB {En' + m] 



dy e'^yiniy) = V27r 



(4.20) 



(4.21) 



In' 



qy_ 



n 2 



-2V2n'eBg^In' 



+ 

In'-l 



In' 



qy_ 

feB 



/eB 



(4.22) 



Thus, 



<9 



V 



- g\)mnmpk ■ A - 2gvgA{P ■ A P' ■ k - P' ■ A P ■ k) 



+s(2gvgAmn{P' ■ Ak^ + P' ■ kA/) - A^k ■ R + k^A- R 



(4.23) 



We now choose the axes such that the 3-momentum of the incoming neutrino is in the x-z 
plane. We will also assume that |P'| <^ nip for the range of energies of interest to us. In that 
case, it is easy to see that the terms involving V^n'eB drop out, and we obtain 



2 2tt mnlTlp 



X — 



eB E„' + m 



{9v + 35i)wAo + (5y - 9\)kzAz 
+2gAS[{gv - gA)t^^z + {gv + gA)kzAi 



(4.24) 



We now put this expression into Eq. ()4.13|) and calculate the total cross section by performing 
the integrations over different final state momenta appearing in that formula. First we integrate 
over and P'^. These appear only in the momentum conserving 5-function. Integration over 
them therefore just gets rid of the corresponding 5- functions. For the integration over p'^, we 
refer to Eq. ()2.12|) . Since the center of the oscillator has to lie between —\Ly and \Ly, we 
conclude that —\LyeB < p'^ < \LyeB. Thus the integration over p'^ gives a factor LycB. 
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Putting back into Eq. H4.13() and using V = L^LyLz, we obtain 
G} S{Q + uj-En') 



da 



An 



{gv + 35a)wAo + [gv - 9A)kz^z 
+2gAS[{gv - gA)uj^z + {gv + gA)kzM 



dP'dp'z, (4.25) 



where Q is the neutron-proton mass difference, rUn — rup. 

We next perform the integration over Py. In the integrand, it occurs only as the argument 
of the functions /„ and In-i- The functions /„ are orthogonal in the sense that 



da In{a)In'{a) 



eB 6nn' 



(4.26) 



This property can be used to perform the integration over Py. We have already remarked that 
the term proportional to \/2n'eB in Eq. (|4.22() does not contribute. From other two terms, we 
obtain 



eB 
eB 



where 



(pf -pf)(i-v,o)+(pf +pjn 

gn' = 2 — 5n',( 



(4.27) 



(4.28) 



gives the degeneracy of the Landau level. Notice the appearance of the Kronecker delta, 5n',0) in 
the expression of Eq. ()4.27|) . The reason for this is that, while two terms of Eq. ()4.22|) contribute 
in the integral for n' ^ 0, only one of them contributes for n' = since /_i = 0. 

The final integration is over p'^ . Writing the argument of the remaining 5- function in terms 
of p^,, we find that the zeros occur when 



Therefore, 



6{Q + u;- En' 



p'^=p'^ = ±^J{Q + Lo f - m2 - 2n'eB . 
Q + uj 



y/{Q + uj)'^ -m'^ - 2n'eB 



= {6{p',-p'^)+6{p'z-p'_)). 



(4.29) 



(4.30) 



In the integration, the terms proportional to p'^ in the integrand receive equal and opposite 
contributions from the two 5 functions and cancel. For the other terms, independent of p!,, both 
the contributions are equal. So we obtain 



an' 



eBG} 
2tt 



^i'[{gv + ^9 a) + 2gAS{gv + gA) cos 6^ 
+Sn',o[{gv - 9a)cos9 + 2gAS{gv - 9a)] 



Q + uj 



(4.31) 



y/{Q + uj)'^ -m^ - 2n'eB 
where we have defined the direction of the incoming neutrino by the angle 9, with 

k^=uj cos 9. (4.32) 
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In Eq. ()4.31|) . we have denoted the cross section by (t„/ because the electron ends up in a 
specific Landau level m! . The total cross section is then given as a sum over all possible values 
of n', i.e., 



(y = ^ cTn' = ^ ^ X! 9n'{{gv + 35!) + '^9AS{gv + Qa) cos 

n'=0 ^ n'=0 

+^n'fi{{gv - 9 a) COS e + 2gAS{gv - gA)] 



^ + " (4.33) 



y/{Q + w)^ - - 2n'eB 

The possible allowed Landau level has a maximum, n'^^^, which is given by the fact that the 
quantity under the square root sign in the denominator of Eq. (|4.33|1 must be non-negative, i.e., 

<ax = int I ^ [(Q + ^)2 - m^] I . (4.34) 

Eq. (|4.33|) gives our result for the cross section of the inverse beta decay process. Some 
properties of this formula are worth noting. 

For unpolarized neutrons, S = 0, the cross section for n' 7^ does not depend on the direction 
of the incoming neutrino. The same is not true if the electron ends up in the lowest Landau 
level. The cross section will be asymmetric in this case. 

All terms in the cross section which depend on S have a common factor g^. The reason 
is that, if gA were equal to zero, the interaction in the hadronic sector would have been spin- 
independent. 

If the final electron is in the lowest Landau level and the initial neutrino momentum is 
antiparallel to the magnetic field, Eq. ()4.33() shows that 



2 r 



%i(i-5) 



^ + " (4.35) 
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Note that the vector coupling gy does not contribute to the cross section in this limit. This can 
be understood easily. The neutrino spin is along the +z direction whereas the electron spin in 
the lowest Landau level must be in the —z direction. Thus there is a spin-flip in the leptonic 
sector. Conservation of angular momentum then implies that there must be a spin-flip in the 
hadronic sector as well. In the non-relativistic limit for hadrons that we have employed, this 
can occur only through the axial coupling. 

If further we consider totally polarized neutrons, i.e., S = 1, we see that ctq vanishes. Again, 
this is a direct consequence of angular momentum conservation. Since both initial particles have 
spin up, angular momentum conservation requires both flnal particles in spin up states as well. 
But the spin-up state is not available for the electron in the lowest Landau level. 

It is instructive to check that the result obtained in Eq. (|4.33|) reduces to the known result 
for the fleld-free case. The contribution specific to the zeroth Landau level vanishes in the limit 
B ^ owing to the overall factor of eB. The other terms also have the factor eB, but in this 
case we also need to sum over infinitely many states. This gives 

cBG^a [or, 

a = {gy + Sg^) + 2gAS{gv + Qa) cos 9 

vr L 




Q+uj Q + , 



(4.36) 
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Figure 1: Enhancement of cross section in a magnetic field for an initial neutrino energy of 10 MeV. Different 
panels show the results for different net polarizations of the neutrons. The solid and the dashed lines correspond 
to the initial neutrino momentum parallel and antiparallel to the magnetic field. Each curve has been normalized 
by the cross section in field-free case for the same values of S and cos 9. The horizontal dashed lines represent 
unity in the vertical scale. 



For — > 0, the last term vanishes, and we can identify n'^^^^ as the integer for which the 
denominator of the summand vanishes. Thus we obtain 



eBG' 



a 



IT 



(gv + ^9a) + '^9AS{gv +gA)cos0 / 

J JO 



dn' ^ + " 



y/{Q + Lo)'^ -m^ - 2n'eB 



TT 



(gt^ + 3gi) + 2gAS{gv+g a) cose (Q + u;)J {Q + cof - m\ (4.37) 



which is the correct result in the field-free case. 

In Fig. ^ we have plotted the ratio of the cross section to its corresponding value at i? = 
as a function of the magnetic field. The plots have been done for unpolarized {S = 0) as well as 
totally polarized neutrons along {S = 1) and opposite {S = —1) to the magnetic field, with the 
initial neutrino momentum parallel and antiparallel to the magnetic field. For S" = — 1, we find 
that neutrinos parallel to the magnetic field have a smaller cross section than those antiparallel 
to the field, and the difference is pronounced for large fields. For S = 0, the situation is just 
reversed. For S = 1, if the magnetic field is high enough so that n'j^^^ = 0, we see that the 
cross section vanishes for neutrino momentum antiparallel to the field. The reason for this has 
already been discussed. 
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Figure 2: Cross section for unpolarized neutrons as a function of the magnetic field for a flat energy dis- 
tribution, normalized to the cross section in the field-free case. The initial neutrino momentum is along the 
magnetic field, and energy is 10 MeV. The energy spread uj2 ~ is 0.2 MeV for the left panel and 2 MeV for 
the right panel. 



5 Consequences of neutrino energy spread 

The enhancement factor in Fig. ^ shows some spikes. They appear at values of the magnetic 
field for which the denominator of Eq. (|4.33|) vanishes for some n' . For field values larger than 
this, that particular Landau level does not contribute to the cross section. To the right of the 
final spike that appears in the figure, only the zeroth Landau level contributes. In other words, 
the final electron can go only to the lowest Landau level for such high values of the magnetic 
field. The exact value of B for which this occurs depends of course on the energy of the initial 
neutrino. 

We need to make an important point about these spikes. Each spike in fact go all the way 
up to infinity. The finite height of a spike in the figure is an artifact of the finite step size taken 
in plotting it. 

In reality, of course, a cross section cannot be infinite. In the present case, this is related to the 
fact that the initial neutrinos cannot be exactly monochromatic due to the uncertainty relation. 
There must be a spread in energy, which can be represented by a probability distribution ^(f-i^), 
defined by 

j duj $(cj) = 1. (5.1) 

In that case, the cross section in a real experiment should be written in the form 

a = j doj $(u;)cj(cj) , (5.2) 

where (j{u)) is the expression derived in Eq. (|4.33j) for a single value of energy. 

As an illustration, we consider the case of unpolarized neutrons {S = 0), and take a flat 
probability distribution of initial neutrino energy, viz., 

^ if CJl < W < W2) 



otherwise. 
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Then the integration of Eq. ()5.2() gives 



a 



eBG} 



F{uj2) - F{uji) 



where 



n'=0 



9n'{9- 



V + '^9 a) + ^n' 



i9v - 9 a) cos ( 



m 



2 - 2n'eB . 



(5.4) 



(5.5) 



with n'j^.^^ determined by Eq. ()4.34() . Fig.|2shows the variation of this quantity with the magnetic 
field for cos 9 = 1. In this figure, we normahze the cross section by B 
energy distribution of Eq. (|5.,3j) , which is 



cross section with the 



a{B = 0) 



+ 3gl) / 



(Q + W2) 



m 



3/2 



(Q + ^l)' 



m 



3/2 



(5.6) 



37r{uj2 — i^i) V 

Keeping the central value of neutrino energy as 10 MeV as before, we have drawn these plots 
for two different values of the spread, as mentioned in the caption. For the smaller value of the 
spread in particular, the graph looks very similar to that drawn in Fig. ^ but the difference 
is that now the height of the spikes denote the actual enhancement, and is not an artifact of 
the plotting procedure. For the higher value of the energy spread, we see that the spikes have 
smoothened out. 



6 Comments 

The calculation of the cross section for inverse beta decay process has been performed earlier 
by several authors [HI El- They assumed that the matrix element remains unaffected by the 
magnetic field, only the modified phase space integral makes the difference in the cross section. 
The results they obtained is the same as the term proportional to gy + 3g\ that we obtained. 

We have used exact spinor solutions in a uniform magnetic field to calculate the cross section 
of the inverse beta decay process for arbitrary neutron polarization. We find that, even for 
unpolarized neutrons, the cross section depends on the direction of the neutrino momentum. 
This asymmetry is not surprising since the background magnetic field makes it an anisotropic 
problem. A similar asymmetry has been noted for URCA processes where the neutrino is in 
the final state. The anisotropy for the S = case comes only from the n' = contribution 
However, for n' ^ 0, there is a cancellation between the two possible states in a Landau level 
which washes out all angular dependence in these levels, provided the neutrons are unpolarized. 
The asymmetry in cross section will therefore come only from the n' = state and its amount 
will depend on the relative contribution of this state to the total cross section. If the magnetic 
field is so high that only the n' = state can be obtained for the electron, the asymmetry will 
be large, about 18%. For smaller and smaller magnetic fields, the asymmetry decreases with 
new Landau levels contributing. For polarized neutrons, however, there is an asymmetry even 
in the field-free case. In presence of a magnetic field, the asymmetry will in general depend on 
the magnetic field, as it appears from various plots in Fig. ^ 

This fact can have far reaching consequences for neutrino emission from a proto-neutron 
star. It has been discussed in the literature that the presence of asymmetric magnetic fields can 
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cause asymmetric neutrino emission from a proto-neutron star However, our calculations 

show that even with a uniform magnetic field, neutrino emission would be asymmetric because 
of the cos ^-dependent terms in the cross section. 
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